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SUMMARY 


The method described is an inverse one; the shock shape is chosen 
and the solution proceeds downstream to a body. Bodies blunter than 
circular cylinders are readily accessible^ and any adiabatic index can 
be chosen. The lower limit to the free-stream Mach number available in 
any case is determined by the extent of the subsonic fields which in 
t-orn depends upon the body shape. Some discussion of the stability of 
the numerical processes is given. A set of solutions for flows about 
circular cylinders at several Mach numbers and several values of the 
adiabatic index is included. 


INTRODUCTION 


A simple and fast method for calculation of the flow field between 
a blunt body at zero incidence to a supersonic stream of ideal gas and 
its attendant detached shock wave was presented in reference 1. Further 
work^ including a catalog of solutions for the case of bodies with axial 
symmetry^ was reported in reference 2- These solutions were primarily 
applicable to the subsonic part of the flow field extending from the 
line of symmetry^ and ending in a region where the flow was supersonic . 

The sonic line had no particular significance in this work. 

The problem of supersonic flow around a blunt nose has been one of 
considerable interest in aerodynamics^ both from the standpoint of 
analytic approximation techniques (see refs. 3 and h-) and from that of 
numerical solutions of the full inviscid equations of motion. A survey 
of the then existing literature can be found in reference Ij reference 4 
should also be seen. Current interest in this problem includes the desire 
to determine a solution for a complete configuration. This entails a 
continuation of the results found for the subsonic and transonic regions 
at the nose of the body around the shoulder and on into the purely 
supersonic flow downstream. Several methods for acconrplishing this have 
been proposed in the literature (e.g.^ refs. 5^ 7)^ and some of these 

even include consideration of equilibrium and nonequilibrium real-gas 
effects. The importance of the results of such calculations is so great 
however^ that it was thought worthwhile to develop still another technique 
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for the solution. The present work is therefore a step toward the 
eventual calculation of flow about a whole configuration. In it^ certain 
devices whose desirability was indicated by examination of the results 
of references 1 and 2 are used. Like the method of reference 1^ it is 
indirect in that the calculation begins at a shock wave and proceeds 
inward until a body shape is determined. If this body is not of the 
required shape ^ a different shock is used^ and the process repeated until 
the desired body shape is attained. 

One difficulty that arose in the work of reference 1 had to do with 
the particular coordinate system chosen there. Specifically^ this 
trouble arose when an attempt was made to calculate the flow behind a 
shock wave whose section was an oblate ellipse (shock for which Bs > 1; 
see ref. l) . Thus^ it was decided to abandon the otherwise convenient 
set of orthogonal coordinates devised by Van Dyke (ref. l) and to use an 
oblique set in which such a difficulty could not arise. 

Another point of interest is in connection with the inherent insta- 
bility to random error that is present in attempts to determine numerically 
the solution to a partial differential equation of elliptic type starting 
from initial^ or Cauchy type^ data. This is an ill-set problem of the 
type discussed by Hadamard in reference 8. An example of the effect of 
this instability is given in reference 1. Therefore^ one must compromise 
between truncation error (or the error introduced when a difference analog 
is used for differentiation) and round-off error by using as large a mesh 
size as feasible. This resuilts in some sparsity of data in the upper 
transonic region where a characteristics solution could be started. A 
smoothing operation was introduced into the calculation in an effort to 
get around the instability difficulty. 

Finally; an attempt was made to increase the capability of the method 
to produce desired body shapes by incorporating a second parameter into 
the specification of shock-wave shape. 

Only two-dimensional results are discussed herein^ but the machine 
program is set up so that axially symmetric results are available by 
changing one n-umber in the input data. Some checks with the results of 
reference 2 were made and appeared satisfactory. It should be noted that 
the present results carry the solution only slightly past the sonic line; 
and do not continue into the supersonic region downstream. 

SYMBOLS 


A parameter affecting shock-wave shape (eq. (2a)) 

tluntness of shock or body; respectively (see eq. (2b)) 


ro U) > 


ro u) > 


3 




Cp^Cy constant specific heats at constant pressure and volume^ 
respectively 

f(t) . (eq. (5b)) 

p/ 


P 

Pst 

Rg^Rt) 

S 


s^t 




^^7 

x(y) 


7 

6 

A 

9,0^ 


P 

Ilf 

(jO 


free stream and local Msich numbers ^ respectively 
pressure 

p u ^ 

*^00 00 

value of p on body surface 
value of p at stagnation point 

radius of curvature at nose of shock and hody, respectively 
entropy 

oblique curvilinear coordinates (eq. (l)) 

X and y components of velocity, referred to free-stream 
speed Uoo 

velocity of free stream, along x axis 
Cartesian coordinates, referred to Rs 

functional representation of curve which represents a shock 
wave (eq. (2a)) 

ratio of specific heats, ^ 

Stream angle, 6 = tan'^ 1 

shock Standoff distance 

inclination of tangent to shock or body, respectively (see 
sketch ( c ) ) 

density, referred to free-stream density p 

stream fimction, referred to pcoUoo 

1 

y 
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COORDIKATE SYSTEM 


To avoid some nonanalytic 'behavior of the solution vhich was 
associated with the orthogonal curvilinear coordinate system of reference 
1^ a different curvilinear coordinate system is used here. The require- 
ment of orthogonality will he dropped. A typical example of the coordi- 
nate system is shown in sketch (a) . 

Here the curve x = X(y) represents a 
possible shock-wave shape ^ and the new 
coordinates are simply measured off as 
x = x(y) shown from this original curve and from 

! the axis of symmetry. Thus^ set 


s- 0 


As 


3As 5As 



s = X - X(y) X = s + X(t) 
t = y y = t 

where x^y are Cartesian coordinates made 
dimensionless hy division hy the radius 
of curvature Rg of the shock wave at its 
apex. Lines of constant s and t are shown 
in sketch (a). The (s^t) system is not 
orthogonal but the Jacobian of the trans- 
formation is unity so that there is 
nonanalytic it y introduced by the coordi- 
nate transformation (l) . 


( 1 ) 


0 curve x = X(y) which represents 

a shock wave should have a derivative at 

each pointy and its tangent should be 
vertical where it intersects the symmetry 
Sketch (a) or stream axis . This latter requirement 

is made because only flow fields with 
vertical symmetry are to be considered herein. Many possibilities are 
therefore available;, but because of their simplicity^ conic sections are 
chosen as the basic shapes. (Another attractive possibility^ however^ 
is the power law shapes.) The form chosen for X(y) is now 


X = x(y) = (l + Ay2) 


1 -Jl - Bsy2 
Bs 


(2a) 


A 

3 
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Here ^ Bg is the bluntness parameter introduced in reference 1, and A 
is a second shape parameter which may be used to gain more flexibility 
in the choice of shock-wave shape . The type of conic used is determined 
by the bluntness Bgi 


p 
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Bg < 0 
Bg = 0 
1 > Bg > 0 
Bg = 1 

Bg > 1 


hyperbola 

parabola 

prolate ellipse » 
circle 

oblate ellipse 


(a?) 


In cases where Bg > 0 , it is Just the square of the axis ratio of the 
ellipse . 


An idea of the effect of the parameter A upon the shape of the 
curves chosen to represent the shock waves can be gained from sketch (b) . 
The basic c\irve is a prolate ellipse with Bg = O .5 and the other curves 
show the effect of positive and negative values of A. 


EQUATIONS OF MOTION 

The equations governing the steady 
flow of an ideal, nonviscous ^ nonheat- 
conducting gas in two dimensions can be 
written (subscripts denote differentia- 
tion) 


Continuity: 


(pU)^ + (PV)y = 0 


Momentum: 


UUx + VUy + P^^Px == 0 
UV^ + VVy + p“^Py = 0 


Energy: 


uSx "^^y ^ 


(3a) 


(3b) 


(3c) 



Sketch (b) 


Initial conditions for starting the numerical solution of equations ( 3 ) 
are obtained from the oblique-shock relations. First, however, equations 
( 3 ) will be cast in a form appropriate to the coordinate system discussed 
in the previous section. 


The continuity equation (3a) can be eliminated by the introduction 
of a (dimensionless) stream function t(x,y) such that 

pu = ty ; pv = (k-) 

Since the energy equation (3c) implies that the entropy S is constant 
along streamlines, it follows that S is a function of i|/ alone j 
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s 


= CylOg 


^ = CylOg f(t) 


(5a) 


Th.e relation 


~ = f(t) (5b) 

can te used to eliminate pressure p from the momentum ec^uations (3'b) . 
The results of all this^ and of replacing velocity components hy the 
corresponding derivatives of the stream function^ are as follows: 




rp Px by / 

llf 2 _X ^ ^ ^ ^ 


7+1 


f) 


Px , ff 

^ It 


= 0 


(6a) 


X 


■'l^y^xx + Vxy + Vy ~ 


Pv 

♦x= f * (P 


7+1 


f) 


, ff 
P V V f 


= 0 (6b) 


A 


It might he emphasized at this point that though there appear to be more 
unknowns than equations^ in reality the entropy function f is a function 
of the stream function This follows from the energy relation^ equation 

(3c) • The principle of constant entropy along streamlines holds for equil- 
ibrium flows of a real gas as well as of a perfect gas^ so that equations 
(6) apply to this situation also. The determination of the functional form 
of f(^) will be considered later in the section on boundary conditions. 
However^ it can be said that an analytic form for f(^) can be found only 
for the perfect-gas casej for a Feal gas^ numerical tables or finactional 
approximations thereto must suffice . 

The density derivatives can be eliminated from equations (6) to 
give a single second-order equation in the stream function To do 

this^ Crocco^s relation between vorticity and entropy gradient can be 
used^ and the result is 



2uv 




xy 


+ 


V cV ’yy 7(7-1) 


where c is the local sound speed divided by There did not seem 

to be any advantage gained by working with this single equation rather 
than with the pair (6a) and (6b) . 

In the solution^ the body contour will be located as the locus of 
zeros of the stream function. Unfort*unately^ vanishes all along the 
symmetry axis since this is the streamline containing the stagnation 
point on the bo(Jy. To avoid this^ it is convenient to introduce a 
modified stream function w(x:^y)^ defined by 


oo t>- OJ 


7 


U) 

■which, vanishes only on the hody . 
smaller variation than does in the region of interest. 

In terms of equations (6) become 


= £ 
y 

This new function has also a generally 



Px Py 


(u + yWy)j + y“yx) - + yWy) ~ + yw^ ^ 

-y^U^Wyy + (p7+^f) 


*■ 


^ X /f’ 


PA Py 

(w + yi*^)! -y^xx + y^^x + y^^xC^^x + y^xy) - y^%^ — 


(P^^^f) 


7 n + (U + yu^)^... 


0 (8a) 


0 (8b) 


It remains to transform equations (8) according to the variable 
change of equations (l) . The results of this are the equations 


tX'(oj + tmt)wss + J(w + twt)[u + t(wt - X’Wg)] 



Ps 

P 


— tm^g(u) + tm^ + tX'UJg) - t^t^g( " X*’Wg) + C0g[u> - t ( “ X'Wg)] 


+ tWg[ w + t (m^ - X’Wg) ] — + (p^'*'^f )t(*)gr— - 1 

and 

t(u) + twt)0Jss - [tWs(w + ti^) - 7X' (p^'*'^f) ] 

= tWg(Ug + tc^s - tWg + (p^^^f) 4 ^ + [U + t(w^ 


(9a) 



(91^) 


These are the equations whose numerical solution is sought- Before the 
method of solution is considered, the initial conditions will be discussed. 
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INITIAL CONDITIONS 


If a step-by-step, or marching, method of nvuaerlcal solution is used 
in the present problem, there are the possible alternatives of starting 
either at the body or at the shock. The former, or direct, method has much 
to recommend it since the body shape is nearly always prescribed. However, 
this is not enough information with which to start a solution process, so 
it is necessary to assume, for example, a pressure distribution on the body 
surface. The integration then proceeds until a shock wave can be inserted 
leading to a uniform free stream. If this cannot be done, the pressure 
distribution is modified and the process repeated until a shock wave can 
be successfully inserted. On the other hand, by starting from the shock, 
all the necessary initial data are known immediately and the integration 
can proceed until the body is reached. The iterations in this case are 
then made with altered shock shapes to get a desired body. Of the two, 
this indirect method seems somewhat the simpler to execute and is the one 
to be used here. Furthermore, if consideration is to be given to equil- 
ibrium, and, especially, to nonequilibrium flows of a real gas, then the 
development of the flow field from shock wave to body seems the more 
natural one to consider. Of co\irse, the free-boundary aspect of this 
problem could be avoided if the work were done in terms of, say, t/y, or 
w, and y as independent variables. 


The initial conditions on the shock wave s = 0 will now be deter- 
mined. Sketch (c) shows the shock wave and the local shock angle 0, 

determined from 


y 



tan 9 = 


dx 


2 ^ 

X' 


( 10 ) 


or 

sin^9 = 


1 + X 


i2 ^ 


cos^8 = 




1 


The value of the stream function t in 
the free stream ahead of the shock is ~ Jy 
hence u) = 1 in the free stream. Since the 
stream function is continuous at the shocks 
this is also the starting value of w just 
behind the shock; 

w(0,t) = = 1 (lla) 

From this^ it follows that 


Wt(0,t) = = 0 


(lib) 


<1 on IT— cvj 



ro 04 
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In order to find the starting value of 


t 


X 


Wg, use 

pv 

y 


so that the starting values of p^v need to he known. From reference 9 , 

fn') M^sln^ - 1 fn') (/ + l)M^sin^0 

= 2 cot 0 — ; p''°^ = — ^ — 

(7 + l)M^ 2 + (7 - l)M^sin^0 

It follows that 


u)g(0,t) = 


sxn 


-0(M^ 


sxn 


-0 - l) 


- — cot 0 

t 2 + (7 - l)M^sin^0 


(lie) 


Derivatives along the shock are found with the relation 


= A M = _ X"(t) d^ 

dtj d0 dt 1 + X’^(t) d9 

shock 

The initial values of the quantities Pt/P 

found with this formtilaj they are 


U), 


^tt 


(o) 




st 


(o) 


sin 20 M^sin^0 - 1 


1 

+ - 


X' 


2 + (7 - l)M^sin^0 t 1 + X 


2 cos 20 


(lid) 

M^sin^e - 1 

2 + (7 - l)M^sln^0 


(7 + l)M^sln^20 
[2 + (7 - l)M^sin^ 0 ]^ 


(lie) 


X’ 


1 + X 


_ fa eot 0 - llX-ilMliXSL- 

L 2 + (7 - l)M^sin^0- 


(llf) 


The only remaining quantity whose initial value needs to he known 
is the entropy function f(t). The functional form of f must also he 
found. At the shocks one finds from^ for example > reference 9 that 
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’(o)(^) = 


(o) 


27M^sin20 - (7 - l) 

7(7 + 1)M^ 


2 + (7 - l)M^sin^e 
L (7 + l)M 2 sin 20 


2an7 


(l2a) 


27M^ - (7 - l)(l + X»^) 
7(7 + + X'^) 


'(7 - l)M^ + 2(1 + X’^)' 


(7 + l)M^ 


(12b) 


Now as to the functional form of f one reasons as follows . At the 
shocks f is known as a function of y (or t) through the slope function 
X'(t) as in equation (l2b) . Also^ the stream function hence w, is 
known as a function of t along the shock. Hence ^ along this line^ the 
functional form of f(t) is known^ and this functional form should persist 
throughout the flow field. Ihis now allows the calculation of the initial 
values of the derivative f'/f. Since 

_d ^ dX'^ d 
dt dt dX'^ 

by logarithmic differentiation one has 


9“' 


dX'^(\|f) 

di|f 


7-1 


1 + X'^ 27M^ - (7 - l)(l + X'^) 


2(1 + X'^) + (7 - l)M^J 


(12c) 


The derivative dX'^/dt can be evaluated when the shape of the shock 
wave is prescribed. For the family of curves of equation (2a) 


^ /'2AX(y) , 1 + Ay2 

x'(y) = y + 


X"(y) = 


2 My) 




+ Ay2 M - Bgy^/ 

1 + 5Ay2 , „ 1 + Ay2 


1 + Ay2 s/1 - Bsy2 


+ Bsy2 


(1 - Bsy2)3/2 


(13a) 


(13b) 


There is now available enough information to start the numerical 
solution of equations (9) • 


NUMERICAL PROCEDURE 


The differential equations to be integrated (eqs- (9)) are seen to 
be in the form 




+ B 2 


P 


Cl 

Ca 


(l^a) 

(l4b) 


<; 00 oj 
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where the coefficients A±^ Ci involve only quantities whose values 
are known at the shock wave;, or initial line^ s = 0. Then^ considering 
equations (l^) as linear algebraic equations for Wgg and Pg/p^ these 
quantities can readily he determined at s = 0. There is then sufficient 
information to write truncated Taylor expansions for w and p in terms 
of the independent variable s^ the coefficients of which are known as 
functions of t* These expansions are used to project w and p to the 
next coordinate line s = and the set of values of co and p can 
then be differentiated numerically with respect to y. The whole process 
can be repeated now^ expanding from s = to s = and so forth- 

Sketch (d) shows the arrangement of points used in the computations . 
Notice that^ to avoid difficulty at t = 0^ 
the grid is displaced At/2 off the symmetry 
axis. Results can later be extrapolated to 
y = 0^ making use of the vertical syminetry 
of the flow field. The present conputing 
procedure is set up to allow for a maximum 
of 50 points to be chosen along the initial 
line s = 0 . The interval At and number 
of points (< 50 ) chosen to start the process 
is arbitrary^ except for the following con- 
siderations . The conic section^ or other 
analytic curve ^ which is chosen to represent 
a shock wave will;, in general^ fail to 
represent it well at some distance from 
the symmetry axis. Thus^ if the assumed 
shock is taken to be a circular arc^ then 
only so much of the circle as has a slope 
tan 0 (see sketch (c)) which is 
> (M^ “ 1 )”^'^^ can be used; that is^ 
the curve representing a shock wave must 
be inclined to the synametry axis at an 

angle greater than the free-stream Mach angle ^ otherwise the shock 
relations used to calculate initial conditions do not apply. This 
situation is automatically adjusted in the machine program; the machine 
will drop sufficient points to ensure that the condition above is met . 

The above procedure can be repeated until at some step^ say the 
(n + 1 )^ one or more negative values of ^ appear along the coordinate 
line s = (n + l)As. Since it is clear that the zero stream- 

line for certain values of t (or y) lies between the lines s = nAs 
and s = (n + l)Z^. The location of this zero is then determined from 
the known data along s = nAs by solving for 6s from the equation 
expressing that the truncated series for the stream function vanishes: 

w(s + 6s^t) = 0 = w(nAs^t) -h 6st^s(ixAs^t) + ^ 6s^Wgg(nAs ^t) 

At each value of t^ the body is now located at s = nAs + 6s ; and the 
(s^t) coordinates are easily put in terms of the (x^y) system by means 


y 



Sketch (d) 
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of equations (l) . It can be seen from 
sketch (e) (where a sample of the variation 
of ^ with s at fixed t^ say^ 
t = [m - ( 1 / 2 ) ]At^ is shown) that the 
numerically greater of the two roots of 
the last equation is the proper one to 
choose here . After the ^ 0 streamline 
has been founds the extrapolation process 
begins again and the various quantities 
^s^ ^t^ etc.^ are determined at the 
coordinate line s = (n + 2 )^ . It will 
be noted by the machine that more values 
of ^ have changed sign in this step^ and 
Sketch (e) a new Ca 3 = 0 streamline will be confuted. 

This alternation of extrapolation and body 
determination now continues until the machine is stopped by some form of 
misbehavior of the solution. The usual form of misbehavior is found in 
the excessive growth of the second derivatives of which occurs because 
of the penetration of the body surface after the first body calculation. 
Since there must be some sort of sourcelike singularity inside the zero 
streamline^ calculation within this region will soon fail because of 
excessive growth of the stream function and especially its derivatives. 

Sketch (f) illustrates some of the 
quantities involved in the body deter- 
mination process. 

In solving the linear equations 
(l4) for and Pg/p^ it can happen 

that 


AjlB^ *“ “ 0 


in which case the solution will fail. 

To explain this^ one can derive the 
characteristic form for equations (l4) 
and show that the above condition indicates 
s that at the point in question^ the coor- 
dinate line s = constant is tangent 
to a characteristic line and the failure 
Sketch (f) of the extrapolation procedure is justified. 

This can therefore only occur in regions 
of supersonic flow^ and can generally be avoided by a slight change in 
the interval size or At. 

Calculation of the Plow Variables 

The density p is calculated at every point of the mesh. As in 
reference 1^ an external iteration is performed on density to improve its 
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accviracy. Since the entropy function is also known at each mesh pointy 
calculation of the pressure is immediate^ as 

p = p^f(t) 

When the pressure at each point is known^ the local Mach number is found 
hy means of the formula (see ref. 9) 




2 ^'st 
7 - 1 V Pst 



2 r 2 + (y - l)M ^ p 

7 - 1 L 27M^ P 



(15) 


For cases in which the adiabatic index 
(15) is modified to read 


7 is taken as 1, formula 


= M^(l - sin^e + 


M'^sin^a 


+ 2 log 


sin^a 


(16) 


where sin^a is the same function of t as in the calculation of f (t) • 

The velocity components in the (x^y) system, u and v, as well as 
the stream angle & = tan~^(v/u) are printed out in the machine calcula- 
tion. Ilhese would "be useful for plotting streamlines^ as well as for 
adjoining a characteristics solution to the present one. 

It was mentioned previously that flow variables are not calculated 
on the symmetry axis^ but start At/ 2 off it* To extrapolate these values 
to the axis^ the first two calculated points and the symmetry in y = 0 
are used to determine a cubic parabola* The formula resulting is 


A(s,0) 




(17) 


where A(s^t) represents the physical quantity in question. 


Characterization of Body Shape 


The question of how to characterize a set of body points resulting 
from a calculation naturally arises. One way would be simply to plot 
these points and conpare the results with a plot of the desired shape . 
This is a necessary procedure if the desired shape is not represented by 





a simple formula. However many nose shapes of interest are conic 
sections; in particular the circle is of considerable Importance as a 
standard solution for comparison purposes. Thus^ it is possible to 
gain some idea about the closeness of the calculated body shape to the 
desired one by a nmerical process . The process used here is the same 
as that of reference 2; if the body is a conic,, its equation is 

y2 = 2 — (x - a) - B^(x - A) (l8a) 


where is nose radius,, the blrntness^ and A the standoff 

distance of the body. The quantity 


2(x - A) R, 


- _ 5k (x - A) 


(I8b) 


2(x-A) 


can be plotted versus (x - A)^ and the points approximated to by a 
straight line (see sketch (g)) . The vertical intercept of the line is 

Rjj/Rg, and its slope is -(l/2)Bb. 

^2 The machine program was set up to 

include this calculation hy a least- 
squares fit^ hut this is not completely 
satisfactory because of occasional 
wildness of points^ particularly near 
X - A = 0. Thus^ one can use the 
machine-calculated values of 
y=/2(x - a) to make one^s own fit. 

In addition^ the program calculates 
y at each body value of x based upon 
the least- square values of and . 
These can be checked with the actual 
values and the closeness of the conic 
appr oximat i on e valuat e d . 



X- A 


Sketch (g) 


<; on CM 
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Sketch (h) 

Sketch (h) shows a generalized flow chart for the calculating 
procedure outlined above • 


INSTABILITY AND SMOOTHING 


An attempt was made to examine the effect of the inherent instability 
with respect to random error which is known to exist when a partial differ- 
ential equation of elliptic type is solved numerically with boimdary 
conditions of the Cauchy type - There is of course always a random error 
in ro-unding off to any number of figures which may be kept in a calcula- 
tion^ but in order to get a clearer notion of this effect^ the starting 
values of w (only) were deliberately perturbed from their correct values^ 
and a solution was determined from these new initial conditions* 

The case taken for this purpose was that which resulted in a circular 
cylinder at M = 4^ 7 = 1*4. The shock bluntness Bs = -0.06. Several 
runs were made with errors of different magnitude introduced^ and each 
error was run at two interval sizes in the s direction. In the cases 
with the larger interval size^ As = 0.032^ six steps were required between 
the shock and the body* This number of steps was considered about optimum 
by Van Dyke in references 1 and 2. The smaller interval size^ As = 0.0l6^ 
required 13 steps from shock to body. 
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Values of the modified stream function ^(s,t) at t = (l/2)At for 
the cases with o)[0^ (l/2)At ] = 1-0000000 and wlo^ (l/2)At ] = 1.0004883 
are shown in figure 1. It is seen that the artificial error has a 
considerable effect on the variation of causing a change in standoff 
distance of about 13 percent when the interval size is = 0.032. The 
differences in the case where As = 0.0l6 are much greater^ and the 
solution was^ in fact^ halted by the error accumulation at s = 0.l60. 
The difference between the solutions with no artificial error but with 
As = 0.032 and 0.0l6 are definitely noticeable; this is presumably due 
to the growth of error caused by rounding to eight figures^ as well as 
to the different truncation errors. The standoff distance differs here 
by something less than 2 percent. Also shown in figure l(a) is the 
variation of oo(s^O.ll). At this vertical distance from the induced 
error;, its effect was not visible for as much of the solution as was 
obtained (As = O.OI6) . 

Pigure l(b) shows the variation of co with t along s = constant 
lines for the case As = O.OI6. Here one can see the sudden growth of 
error after several steps have been taken. The broken lines and circles 
represent the values determined with w(O^O.Ol) = 1.0004883^ and the 
solid lines are the results without deliberate error. The differences 
between the data in the two cases are generally too small to appear on 
the plot until s = 0.112. From this line onward^ the growth is very 
fast. As mentioned previously^ q_uantities on the line s = O.176 could 
not be determined because of overflow in the calculator . 

The following table of values of the modified stream function w 
and of the second derivative will give a more precise pictiore of 

the various effects. The numbers are all roimded to four significant 
figures for ease of perusal^ except for the initial values of oj where 
the errors would not then be apparent. 



w[s,(l/ 2 )At] 

s 

As = 0.032 

As = 0.016 1 

0 

1 .0000000 

1 .0000010 

1.0004883 

1 .0000000 

1.0000010 

1.0004883 

.032 

.8846 

.8846 

.8851 

.8845 

.8845 

.8849 

.064 

.7657 

.7658 

.7660 

.7650 

.7650 

.7647 

.096 

.6407 

.6407 

.639^ 

.6384 

.6384 

.6330 

.128 

.5053 

•5053 

.4977 

.5004 

.5004 

.4660 

.160 

• 35*^1 

.3540 

.3220 

.3446 

.3443 

.1565 

.192 

.1788 

.1786 

.05742 

.1617 

.1600 

— 

.224 

- -03299 

- -03379 

- -5173 

- .06576 

- .07005 

— 

j ^SS 

0.032 

-4.432 

- 4.433 

-4.885 

-4.398 

- 4.400 

-5.359 

.064 

-7.732 

-7.737 

-9.872 

-7.687 

-7.697 

-12.78 

.096 

-12.35 

-12.36 

- 19.48 

-12.37 

- 12.42 

-39.82 

.128 

-18.77 

-18.83 

-47.10 

-19 -02 

-19.28 

-162.8 

.160 

-28.04 

-28.21 

-126.4 

-28.93 

-30.15 

-785.3 

.192 

-43.41 

- 44.04 

-479.4 

-47.08 

-50.67 

— 

.224 

-24.90 

-31.15 

-1113 

-96.54 

-131 .1 

— 




Effect of Error at Shock on P; Ps 


^ = 0.032 

s 

w(0,0.0l) =1.0000000 

co(o,o.oi) =1.0000010 

w(0,0.0l) = 1.0004883 

P 

Ps 

P 

P 

Ps 

P 

P 

Ps 

P 

0.032 

.064 

.096 

.128 

.160 

.192 

.224 

4.6820 
4 . 7734 
4.8498 
4.9134 
4.9641 

4.9987 
5 .0080 

3.114 

2.597 

2.179 

1.796 

1.372 

.7873 

-.2009 

0.8539 

.8774 

.8971 

•9136 

.9269 

•9359 

.9384 

4.6820 

4 . 773 ^ 

4.8498 

4.9134 

4.9641 

4.9986 

5.0079 

3.114 

2.597 
2.179 
1.796 
1.372 
.7863 
- .2087 

0.8539 

.8774 

.8971 

.9136 

.9269 

.9359 

■9383 

4.6821 

4.7737 
4.8506 
4 .9152 
4.9661 
4.9835 
4.7444 

3.117 

2.606 

2.202 

1.835 

1.349 

- .2640 
- 14.68 

0.8540 

.8775 

.8973 

. 914 i 

.9274 

.9320 

.8699 

2ss = 0.016 

• 032 
.064 
.096 
.128 
.160 
.192 
.224 

4.6819 

4.7738 

4.8512 

4.9160 

4.9678 

5.0027 

5.0095 

3-133 

2.628 

2.217 

1.833 

1.392 

.7509 

-.4437 

.8539 

.8775 

.8975 

•9143 

.9279 

.9370 

.9388 

4.6819 
4.7738 
4.8512 
4 .9160 
4 .9679 
5.0027 
5.0086 

3-133 

2.628 

2.217 

1.833 

1.392 

.7433 

- .1^88 

•8539 

.8775 

.8975 

.9143 

.9279 

.9370 

•9385 

4 . 6820 
4 . 7743 
4.8532 
4.9224 
4.9633 

3.138 
2.650 
2.297 
2.018 
- .1636 

.8540 

.8776 

.8980 

.9160 

.9267 


The enormous grovth of error is quite apparent when the second 
derivative ^ss is inspected. It is to be noted that halving the 
interval size has an effect on the original solution that is numerically 
greater than that of tbe artificial error of 1X10"“^ and less than that 
of the artificial error of 5X10"’^. A tabulation of the values of density 
p^ the derivative pg^ and pressixre p^ is also shown. Here the error is 
much less striking^ probably due to the external iteration performed on 
p. From this exanple^ it would appear worthwhile to try iteration on 
the stream function also in future work. 

In an attempt to offset the inherent instability discussed above ^ a 
smoothing process was inserted in the machine program. For simplicity^ 
an available SHAEE routine was used.^ This process causes a recalculation 
of the numerical values of the flow quantities along each s = constant 
line by means of a least- square fit of a third-degree polynomial to seven 
consecutive points. Wild points are discarded and replaced by calculated 
values^ and numerical differentiation based on a three-point formula is 
then performed. This routine is actually somewhat more elaborate than 
is necessary for the present purposes^ and a heavy penalty is paid in 
computing time as a result. When the program runs without the smoothing 
procedure^ seven-point numerical differentiation formulas are used for 

the derivatives with respect to t. In this latter ins tance^ a single 

ifhe SHARE program in question is referenced as CL SMD3;j SHARE distri- 
bution no. 331^ entitled ^'Smooth and Differentiate Unequally Spaced Points.” 
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case is completed in about 1 minute on an IBM 704 computer, but with the 
smoothing procedure operative, the average case runs I5 to 20 minutes. 

In order to assess the effect of the smoothing operation inserted 
in the program, the case discussed above was run at several inter-vals. 

It should be mentioned that the unsmoothed program, in which seven-point 
differentiation with respect to t is used, is set to drop the three 
values of the dependent variables corresponding to the three highest 
values of t at each step. This is done to keep central difference 
formulas applicable throughout, because it was believed that uncertainties 
propagated from these top positions were a major cause for error growth. 
Thus, it is not possible to use a very small interval in connection with 
the uismoothed program, for the number of steps to the body is limited 
to 15 by the point-dropping. This limitation does not exist in the pro- 
gram with smoothing. 

A comparison of the values of the modified stream function w as 
calculated using the smoothed and unsmoothed programs is shown in figure 
2 . Figure 2 (a) shows the variation of w[s,(l/ 2 )At] with s for several 
values of As, and for runs made with and without smoothing. In each 
case, the effect of the smoothing is to raise the level of w over the 
entire course of its variation, without affecting the general form of 
the curve. This is due to smaller numerical values of the second 
derivative ^^ss being found in the smoothed case . It does not seem 
possible to say whether this is affecting the solution process in a 
proper way or not, but it does appear that the closeness of the results 
in cases run with and without smoothing and the established tendency of 
the unsmoothed solution toward enlargement of errors make the likelihood 
of an adverse effect by the smoothing process rather small. A table of 
values of w for the various cases is included so that a closer 
comparison can be made . 


s 

00 

Smootlied 

Unsmoothed 

As = 0.032 

As = 0.016 

li 

0 

b 

0 

L^- 

= 0.032 

As = 0.016 

0 

1 .0000 

1.0000 

1.0000 

1 .0000 

1 .0000 

.032 

.8846 

.8845 

.8844 

.8846 

.8845 

.064 

.7658 

.7650 

.7646 

.7658 

.7650 

.096 

.6409 

.6388 

.6375 

.6407 

.6384 

.128 

.5061 

■5015 

.4988 

.505^ 

.5004 

.160 

•3559 

•3^73 

.3^21 

• 35^1 

.3446 

.192 

.1830 

.1681 

•1591 

.1788 

.1617 

.224 

- .02382 

- .o48o4 

- .06240 

- .03299 

- .o48o4 


<; 00 oj 
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In figiare 2(15)^ stagnation pressure Pg^ and stand-off distance 
fi/Rg are plotted versus interval size for several smoothed and unsmoothed 
runs, (it should he mentioned that in the smoothed calculation^ the 
entropy fmction f (''^) is fixed at its known value on the body when pres- 
sures on the body are computed.) The dashed line on the pressure plot 
represents the known exact value of stagnation pressure. Two additional 
smoothed runs were made, with As = O.OO8 and 0.010 (the case with 
As = 0.008 also appears in fig. 2(a)), and the results were plotted in 
figure 2(b) . The plot of stagnation pressures against interval size As 
seems rather inconclusive since sufficiently small intervals for comparison 
are not available in the unsmoothed case. However, the runs with smoothing 
do approach the known correct value in a rather consistent manner. 

In the case of the plot of detachment distance A/Rg versus As, it 
is noted that the smoothed points lie on a line so that linear extrapola- 
tion to the limit of vanishing interval size (in the s direction) is 
possible. On the other hand, the unsmoothed points for A/Rg lie on a 
line until the Interval size As becomes O.OI5, and this is what one 
would expect; that is, the solutions will begin to be less accurate as 
the nimiber of steps taken becomes too large . Unfortunately, smaller 
values of As cannot be used in the unsmoothed program because the body 
streamline would not be reached. The nearly constant difference in level 
between the smoothed and -unsmoothed extrapolation lines follows from the 
properties of these solutions mentioned above in the discussion of figinre 
2(a) . Smaller intervals for the smoothed case were not run because of the 
large amount of machine time used. For As = O.OO8, 26 steps were taken 
to reach the body. 

Several runs were made with the smoothed program in which a deliberate 
error was forced in the initial data. Figure 3(a) shows plots of w(s,t) 
versus t at several values of s proceeding from the shock toward the 
body. The error in initial data was very large, the value of u) being 
taken as 1.1 at s = 0, and t = 0.01. The propagation of the error can 
be followed along the successive s = constant lines in figure 3(a), and 
its spread and intensification are fairly rapid. The error is, of coiirse, 
enormous for a round-off error, but its large size makes the effect easier 
to follow. Figures 3(b) and 3(c) show w as a function of s for t = 0.01 
and t = 0.21, respectively. The effect of errors in w(0,0.0l) of 0.1 
and 0.01 are shown in figure 3 (b), while only the negligible effect of 
the larger error on w(s,0.2l) is shown in figure 3(c). A deliberate 
error introduced at a higher value of t, for example, t = 0.11, has a 
much smaller effect on the solution, as was remarked by Van Dyke in 
reference 1. The values above and below t = 0.11 are apparently unaf- 
fected, and the bump in the data merely propagates across to the body, 
being lightly dan5)ed rather than growing . 

Finally, it is the author's opinion that the use of smoothing in the 
numerical solution does not affect the results adversely, but has the 
effect of allowing the use of a smaller Interval size than is feasible 
with an unsmoothed program. It should be stated, however, that the 
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instability of the set of equations (9) is not as serious as might have 
been expected. This may be because the end points of each s = constant 
line are dropped, to avoid error vhich could arise from the noncentral 
derivative formula. There is also a strong possibility that the nonlin- 
earity of equations (9) has an alleviating effect on the instability. In 
this latter connection it is interesting to note that the present method 
of numerical solution was attempted for an incompressible flow about a 
circular cylinder, where the governing equation is Laplace's. It was 
found that the Instability was disastrous in this case, even when the 
data were subjected to smoothing. Thus, some connection might be made 
between the "strength of the ellipticity" of an equation and its degree 
of instability toward random error when used in connection with boimdary 
conditions of Cauchy type . 


NUMERICAL RESULTS 


As an illustration of the dependence of such geometric quantities 
as B]3 on the shock par'ameter Bg, several cases 

were run off at M = and 7 = 1.4 with -2 < Bg < + 2. The results 
are summarized in figure 4. The bodies resulting from the various values 
of Bg were not in all instances well fit by conic sections. This was 
particularly true for the smaller values of Bg. However, in these cases, 
the stagnation region was reasonably well approximated by the conic section 
so that the results hold at least for this portion of the bodies. 

The additional geometric parameter A, introduced in equation (2a) , 
has the effect of making very blunt bodies available . Figure 5 illustrates 
this for a free stream at M = 10 and 7 = 1.4. The body corresponding 
to A = 1 is actually slightly concave near the stagnation point. The 
A = 0 body is a circular cylinder, and the A = -1 body is nearly para- 
bolic. In the latter case, the sonic point on the body was not reached. 

In fact, sonic speed behind the shock was not attained on the portion 
shown in figure 5* The pressure distributions at the body surfaces 
for the various cases are also shown, plotted against the vertical 
coordinate y. 

The variation with 7 of the geometric parameters of a body corre- 
sponding to a parabolic shock wave with infinite free- stream Mach number 
is shown in figure 6 where the shock wave and the various bodies are 
plotted. For 7 = 1 at M = 00^ the body and shock axe coincident. The 
calculations were done for 7 = I.O5, 1.1, 1.2, 1.4, I.6667, 2, and ?. 

The table insert in figiore 6(a) gives the numerical values of the 
geometric parameters associated with the several bodies. In addition, 
the variations of and Rb/Rg are plotted against (7 - l)/(7 + l) . 

From the results shown in figure 6, one can see that the parameter B^ 
is not enough to specify a shape; the radius R5 must also be given. 


00 OJ 
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Flow Fields About Circular Cylinders 


The main results of the present calculations are contained in table 
where some characteristics of flow about circulaj: cylinders at various 
free-stream Mach numbers and with several values of the adiabatic index 
7 are listed. The Mach mimb ers used are M = 20^ 10^ 8^ 6 , 4^ 3^ 

and 2 .f and for each of these, the parameter 7 takes the values 6/5, 

7/5^ ao,<i 5/3 • Quantities given for each set of conditions are the 
coordinates of shock, body and sonic line, and pressures on the body. 

It did not prove feasible to give results for free-stream Mach 
nimibers less than 2 for these two-dimensional flows . The greater stand- 
off distances and generally large extent (based on body size) of the 
subsonic field in plane flows as opposed to axially symmetric flows make 
determination at low Mach numbers very difficult. If one confined 
attention to bodies rather blimter than a circular cylinder, solutions 
could probably be found down to ^ch numbers of about 1.2, but precise 
results in this range are not presently available. For the M = 2 cases 
with 7 = 7/5 aJid 5/3^ a well-established sonic line was not found. This 
is presumably because of its large extent through the flow field and 
consequent sensitivity to errors made at the largest values of the 
vertical coordinate t . 

Most of the solutions listed in table I were found using the 
smoothing process discussed above. However, some cases were determined 
before this process was incorporated. Spot checks were made to compare 
the results with and without smoothing, and the closeness thereof in the 
cases checked gave confidence in using the unsmoothed results where avail- 
able. In all cases, the experimental rims to determine the approximate 
value of Bs at a given M and 7 were made with the unsmoothed process 
because of the much greater speed of that program. 

The pressure ratios given in table I are the ratios of local pressure 
to . stagnation pressure Pg^ at the apex of the body. The values of p^^ 
used are those given by the numerical solution, and are listed for each 
case, along with the known exact ■values computed from (see ref. 9) 

1 


The calculated values of pressure ratio at M = 00 may be compared with 
values given by the "modified Newtonian" pressure formula, namely 

P p 

- — = sin^db 

^st 




(19) 
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where is the angle made with the stream direction hy a tangent to 

the "body surface. (For a discussion of Newtonian results^ see ref. ^.) 

The comparison (see fig. 7) shows fair agreement between the calculated 
values and those given by equation (l 9 ) • 

The fortuitous agreement between modified Newtonian and actual pres- 
sirre distributions at 7 = 1 *^ is demonstrated^ though in lesser degree 
than in axially symmetric flow. The behavior of the pressure distributions 
for decreasing 7 is somewhat surprising in that for 7 = 5 / 3 ^ 7 / 5 ^ and 
6 / 5 , the several pressure distributions fall first below^ then above the 
modified Newtonian result. This is seemingly the wrong trendy for as 
7 1 ^ the limiting pressure distribution should be the Newton-Busemann 

one^ where the centrifugal correction has been made. The formula for this 
is 

^ = sin2©^ - I cos^a^ 

•t'st 

(see ref. 4)^ and a plot of it is also shown in figure ?♦ In order to 
determine whether the trend is not yet established because 7 = 1.2 is 
still too far from 7 = Ij additional solutions were run for a circular 
cylinder at infinite Mach number and for 7 = 1.1 and 1.05* The results 
from these runs are plotted in figure 7 ^ and one can see that they do 
indeed fall as they should. Thus, it is seen that 7 must be quite near 
to unity in order that the limiting case of Newtonian flow be a fair 
representation of the actual results. 

The results for M = 20 fall very nearly as those for M = and 
results from a version of equation ( 19 ) corrected for finite Mach nirmber 


p 1 + 7M^sin2&i3 

I* St 1 + 7M^ 


(19a) 
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hardly differ from those for M = oo. For the lover Mach numbers^ however 
agreement between calculated and modified Newtonian results becomes 
increasingly poor^ the calculated results being generally the higher so 
that the centrifugal correction would only worsen the situation. 


Unfortunately the data for two-dimensional blunt-body flows are 
not so plentiful as for axis;y7nmetric ones- However^ some numerical 
solutions are available in references 10 and 11^ and reference 12 gives 
a compilation of e^q)erimentair data. Standoff distances from these refer- 
ences are plotted with the present results in figure 8. The results of 
reference 10 and the present ones differ from the experimental data of 
reference 12 by about equal amounts but in opposite directions. The 
asymptotic value given in reference 11 is matched very closely by the 
present result^ calculated for M - 10^. 
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The results for 7 = foimd here also correlate well on the 
basis of the ’’detachment angle” parameters discussed in reference 13- 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Fields Calif. ^ Apr. 11^ I961 
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TABLE I.- GEOMETRY AED PRESSURE DATA FOR CIRCULAR CYLIEDERS 


Bs = 0.1; Rb/Es = 0-5003; ^Rg = O.1888; ^ 


Pressijire 

ratio 


P/Pst 


1.000 

.996 

.979 

.9^9 

.906 

.852 

.786 

.710 

.625 

.529 

.4-17 


Shock 

y 

X 

0 

0 

.03 

.00045 

.07 

.00245 

.11 

.00605 

.15 

.01126 

•19 

.01807 

.23 

.02649 

.27 

.03652 

.31 

-o 48 i 7 

•35 

.06144 

•39 

.07634 

.43 

.09288 


0 


3 

.44 

8 

.05 

12 

.7 

17 

.4 

22 

.3 


= 0.9193 (0.9197) 


Sonic line 


0.0825 0.405 


M = 20 


Bs = 0.16; Rb/Rg = 0.4991; A/Rs = 0.1908; ps^. = 0.92062 (0.92051) 


1.000 

.996 

.979 

.950 

.907 

.852 

.786 

.708 

.618 

.514 

.389 


M = 10 


Bg = 0.1; Rb/Rg = 0.4860; A/r^ = 0.1942; Pg^ = 0-92299 (0. 92298) 


0 

0 

.03 

.00045 

.07 

.002450 

.11 

.006053 

•15 

.01126 

.19 

.01808 

.23 

.02651 

.27 

.03656 

.31 

.04824 

.35 

.06155 

.39 

.07652 

.^3 

.09314 


0 

0 

.1942 

0 

1.000 

.0860 

. 4 i 4 

.015 

.000113 

.1945 

1.8 

.999 

.12 

.434 

.045 

.001013 

.1963 

5.3 

.991 

.16 

• 453 

.075 

.002813 

.2001 

8.9 

•975 

.20 

•459 

.105 

.005514 

.2057 

12.5 

.952 

. 2 k 

• 453 

•135 

.009117 

.2133 

16.1 

.921 

.28 

.417 

.165 

.01362 

.2230 

19.8 

.883 

. .30 

.384 

.195 

.01903 

•2349 

23.6 

.839 

.32 

.339 

.225 

.02534 

.2492 

27-6 

.788 

— 

— 

•255 

.03257 

.2662 

31.6 

• 731 

— 



.285 

.04070 

.2863 

35-9 

.667 

— 

— 

.315 

.04974 

.3099 

4 o .4 

.598 

— 

— 

•345 

.05969 

•3381 

45.2 

.521 

— 

— 

.375 

.07056 

-3722 

50.5 

.436 

1 

— 

.405 

.08235 

.4152 

56.4 

.338 


— 

•435 

.09506 

— 

— 

— 


: 



























TABLE I.- GEOMETRY MB PRESSURE BATA FOR CIRCULAR CYLINDERS - Continued 

(a) 7 = lA - Continued 

1 


Bg = O.lj Rb/^s = 0 . 478 iij a/Rs = 0.1969; Pst = 0.92533 (o. 92 i; 6 i;) 


Shock 

Body 

Pressure 

ratio 

Sonic 

line 

7 

X 

X 

90° - ©b 

P/PBt 

X 

y 

0 

0 

0 .1969 

0 

1.000 

0.0880 

0 . 4 l 8 

.03 

.0005 

.1978 

3.6 

.996 

.12 

.437 

.07 

.0025 

.2020 

8.4 

.978 

.16 

.454 

.11 

.0061 

.2097 

13.3 

.946 

.20 

.462 

.15 

.0113 

.2210 

18.3 

.901 

.24 

.452 

•19 

.0181 

.2362 

23.4 

.843 

.28 

.415 

.23 

.0265 

.2557 

28.7 

.773 

.30 

.366 

.27 

.0365 

.2803 

34.4 

.691 

— 

— 

.31 

.0482 

.3111 

to. 4 

.598 

— 

— 

.35 

. 06 l 4 

.3499 

47.0 

.491 

— 

— 

.39 

.0763 

.4011 

54.6 

.367 

— 

— 

.43 

.0929 

— 

— 

1 

— 

— 


M = 6 


Bs = 0; Rb/Rg = O.A558J a/Rs = 0.2026; p^t = 0 .9283 (O.9289) 


0 

0 

.2026 

0 

1.000 

.094 

.433 

.03 

.0004 

.2036 

3.8 

.996 

.12 

.447 

.07 

.0024 

.2081 

8.8 

.976 

.16 

.461 

.11 

.0060 

.2161 

l 4 .o 

.942 

.20 

.467 

•15 

.0112 

.2279 

19.2 

.893 

.24 

.460 

.19 

.0180 

.2439 

24.6 

.830 

.28 

.427 

.23 

.0264 

.2646 

30.3 

.755 

.30 

.384 

.27 

.0364 

.2906 

36.3 

.668 

— 


.31 

.o 48 o 

.3234 

42.8 

.569 

— 



.35 

.0612 

.3653 

50.2 

.459 

— 

i 

•39 

.0760 

— 

— 

— 

— 

— 

.43 

.0924 

— 

— 

— 

— 

— 


M = A 


Bs = -0.06; R-b/Rs = 0.4151; A/Rg = 0.2182; Pgt = 0.94123 (0.94054) 


.107 

.464 

.12 

.468 

.16 

.476 

.20 

.477 

.24 

.458 

.28 

.409 

.30 

• 356 


0 

0 

.2182 

0 

1.000 

.03 

.0004 

.2193 

4.15 

•995 

.07 

.0024 

.2242 

9-71 

.973 

.11 

.0060 

.2330 

15.4 

.933 

.15 

.0112 

.2462 

21.2 

.877 

.19 

.0180 

.2642 

27.2 

.806 

.23 

.0264 

.2876 

33.6 

.719 

•27 

.0364 

.3178 

to .6 

.618 

.31 

•0480 

.3570 

48.3 

.501 

.35 

.0611 

— 

— 


.39 

.0759 

— 


.. — 

.43 

.0922 

— 

— 

— 

.47 

.1101 

— 

— 

— 











TABLE I.- GEOMETRY AHD PRESSURE DATA FOR CIRCULAR CYLIEDERS - Continued 

(a) 7 = lA - Concluded 


Bg = -0.26; R^/Rg = 0.3602; A/Rg = 0.2362; Pg^. = 0.9577 (0.95722) 


Shock 

y 

X 

0 

0 

.03 

.0004 

.07 

.0024 

.11 

.0060 

.15 

.0112 

.19 

.0180 

.23 

.0264 

.27 

.0363 

.31 

.0478 

•35 

.0608 

.39 

.0753 

.43 

.0914 

•^7 

.1089 

•51 

.1279 


Body 


0.2362 

.2374 

.2i^30 

.2533 

.2688 

.2903 


0 

4.78 

11.2 

17.8 

24.6 

31.8 


.3191 39.7 
.3579 ^-6 


Pressure 

ratio 


- ^ P/Pst 


1.000 

.994 

.966 

.916 

.844 

.753 

.641 

.510 


Sonic line 


0.127 

.16 

.20 

.24 

.28 

.30 


= -1.5; Rb/Rs = 0.2287; A/Rs = O.266I; = 1.0062 (1.0072) 


0 

0 

.03 

.0004 

.07 

.0024 

.11 

.0060 

•15 

.0112 

.19 

.0178 

.23 

•0259 

.27 

.0355 

.31 

.0464 

•35 

.0587 

.39 

.0721 

.43 

.0868 

.47 

.1026 

•51 

.1194 

•55 

.1371 

•59 

.1558 

.63 

.1754 

.67 

.1957 

.71 

.2168 

.75 

.2386 


.2661 

.2680 

.2770 

.2942 

.3222 

.3679 


0 

7.5^ 

17.8 

28.8 
4i.o 
56.2 

























TABLE I.- GEOMETRY MD PRESSURE DATA FOR CIRCULAR CYLIHDERS - Continued 

(b) 7 = 1.2 


1 M = oo 1 

Bs = 0 .12 j Rb/Rs " 

0.62085 A/Rs = 0.12625 Pst = 

0.95557 (0.95550) 1 

Shock 

Body 

Pressure 

ratio 

Sonic 

line 

y 

X 

X 

90° - 9b 

P/Pst 

X 

y 

0 

0 

0 .1262 

0 

1.000 

0 .o 45 

0.300 

.025 

.0003 

.1267 

2.31 

.998 

.08 

.343 

.065 

.0021 

.1296 

6.01 

.988 

.12 

.386 

.105 

.0055 

.1352 

9.73 

.969 

.16 

.422 

.145 

.0105 

.1434 

13.5 

.941 

.20 

.448 

.185 

.0171 

.1545 

17.3 

.906 

.24 

.459 

• 225 

.0254 

.1685 

21.2 

.863 

.28 

.432 

.265 

.0352 

.1856 

25.3 

.813 

— 

— 

.305 

.0466 

.2059 

29.4 

.757 

— 

— 

.345 

.0597 

.2302 

33.8 

.696 

— 

— 

.385 

.0744 

.2588 

38.3 

.629 

— 

— 

.425 

.0908 

.2928 

43.2 

.554 

— 

— 

.465 

.1088 

.3347 

49.3 

.46o 

— 

— 

M ^ 20 ] 

Be = 0.1l4j Ri^/Re = 

0.60495 

A/Rg = 0.12805 Pst 

= 0.95675 (0.95650)] 

0 

0 

.1280 

0 

1.000 

.046 

.303 

.03 

.0005 

.1238 

2.84 

.997 

.08 

•345 

.07 

.0025 

.1321 

6.65 

.986 

.12 

.387 

.11 

.0061 

.1381 

10.5 

.965 

.16 

.422 

•15 

.0113 

.1469 

i4.4 

.936 

.20 

.446 

.19 

.0181 

.1586 

18.3 

.898 

.24 

.458 

.23 

.0265 

•1733 

22.3 

.854 

.28 

.439 

.27 

.0365 

.1914 

26.5 

.802 

— 

— 

.31 

.0^82 

.2132 

30.8 

.745 

— 

— 

.35 

.0615 

.2385 

35A 

.676 

— 

— 

.39 

.0764 

.2690 

40.1 

.607 

— 

— 

.i^3 

.0929 

.3057 

45.3 

.525 

— 

— 

A7 

.nil 

.3501 

51.0 

.433 

— 

— 

M = 10 

Bg = 0.11 

; R^/Rs “ 

0.59045 

A/Rg = 0.13425 Pst = 

0.95948 (0.95949) 


0 

0 

.1342 

0 

1.000 

.049 

.015 

.0001 

.1344 

1.46 

.999 

.08 

.045 

.0010 

.1359 

4.37 

.994 

.12 

.075 

.0028 

.1389 

7-29 

.980 

.16 

.105 

.0055 

•1435 

10.2 

.967 

.20 

*135 

.0091 

•1497 

13.2 

•945 

. 2 h 

.165 

.0136 

•1575 

16.2 

.920 

.28 

.195 

.0190 

.1670 

19.3 

.890 

— 

.225 

.0253 

.1782 

22.4 

.855 

— 

-255 

.0326 

•1913 

25.6 

.815 

— 

.285 

.o 4 o 7 

.2064 

28.9 

.774 

— 

.315 

.0497 

•2237 

32.2 

.727 

— 

.345 

.0597 

•2435 

35-8 

.679 

— 

.375 

.0706 

.2660 

39-4 

.624 

— 

.405 

.0824 

•2919 

43.4 

.567 

— 

.435 

.0951 

.3222 

47.5 

.500 

— 

.465 

.1088 

•3591 

52.0 

.423 

— 










TABLE I.- GEOMETRY AND PRESSURE DATA FOR CIRCULAR CYLINDERS - Continued 

(u) 7 = 1.2 - Continued 


M = 8 


Bg = 0.06j R^Rg = 0. 571^3 J A/Rg = 0.1384 j Pgt = 0.96170 (0.96175) 


Shock 

Body 

Pressure 

ratio 

Sonic 

line 

y 

X 

X 

90° - 9b 

p/Pst 

X 

y 

0 

0 

0 .1384 

0 

1.000 

0.052 

0.323 

.03 

.0005 

.1392 

2.99 

.997 

.08 

.352 

.07 

.0025 

.1427 

7.00 

.985 

.12 

.391 

.11 

.0061 

.1491 

11.0 

.963 

.16 

.422 

.15 

.0113 

.1584 

15.1 

.932 

.20 

.444 

.19 

.0181 

.1708 

19.3 

.893 

.24 

.452 

.23 

.0265 

.1864 

23.6 

.846 

.28 

.431 

.27 

.0365 

.2057 

28.0 

.790 

.30 

.4o4 

.31 

.048i 

.2289 

32.7 

• 726 

— 

— 

.35 

.06i4 

.2568 

37.6 

•659 

— 

— 

.39 

.0762 

.2907 

42.8 

.587 

— 

— 

.i^3 

.0927 

.3322 

48.5 

.503 

... 

— 


M = 6 


Bs = 0.106; R^/Rs = 0.5555 j a/Rs = o.l47lf; Pst = 0. 96670 (0.96664) 


0 

0 

.1474 

0 

1.000 

.057 

.338 

.015 

.0001 

.1475 

1.55 

■999 

.08 

.360 

.045 

.0010 

.1492 

4.65 

.994 

.12 

.394 

.075 

.0028 

.1524 

7.76 

.983 

.16 

.421 

.105 

.0055 

.1573 

10.9 

.964 

.20 

.438 

.135 

.0091 

.1639 

l4.1 

.940 

.24 

.438 

.165 

.0136 

.1722 

17.3 

.910 

.28 

.391 

.195 

.0190 

.1824 

20.5 

.878 

— 

— 

.225 

.0253 

.1945 

23.9 

.841 

— 

— 

•255 

.0326 

.2087 

27.3 

.797 

— 

— 

.285 

.0407 

.2250 

30.9 

.752 

— 

— 

.315 

.0497 

.2443 

34.5 

.699 

— 

— 

.345 

•0597 

.2663 

38.4 

.643 

— 

— 

.375 

.0706 

.2918 

42.5 

.582 

— 

— 

.405 

.0824 

.3218 

46.8 

.510 


1 

.435 

.0951 

— 

— 

— 




M = 4 














TABLE I.- GEOMETRY AM) PRESSURE DATA FOR CIRCULAR CYLIMDERS - Continued 

(u) 7 = 1.2 - Concluded 


0.i3j R-h/Rr = 0.4238; z^Rg = 0.1977; p,+ = 1.0017 (i.ooio) 


Shocic 


Body 


.0004 

.0024 

.0060 

.0112 

.0180 

.0264 

.0364 

.0479 

.0610 

.0757 

.0919 

.1097 


B3 - -1.15; Rb/^£ 


0.1977 

.1988 

.2035 

.2122 

.2251 

.2427 

.2656 

.2952 

.3336 


Pressure 

ratio 


- I p/R£ 


1.000 

.996 

.976 

.940 

.890 

.825 

.746 

.652 

.541 


Sonic 

line 

X 

y 

0.093 

0.433 

.12 

.442 

.16 

.448 

.20 

.442 

.24 

.4i8 

.28 

•341 


M = 2 


= 0.2652; A/Rg = 0.2454; pg^ = 1.0607 (1.0616) 


0 

0 

.2454 

0 

1.000 

.03 

.0004 

.2471 

6.49 

•990 

.07 

.0024 

.2548 

15-3 

.946 

.11 

.0060 

.2692 

24.5 

.866 

.15 

.0112 

.2919 

34.4 

• 748 

.19 

.0179 

.3262 

45.8 

.591 

.23 

.0261 

•3795 

60.1 

.396 

■ 27 

.0357 

— 

— 

— 

• 31 

.0468 

— 

— 

— 

•35 

.0592 

— 

— 

— 

■ 39 

•0730 

— 

— 

— 

•43 

.0880 


— 

— 

■47 

.1042 

— 

— 

— 

•51 

.1216 

— 



— 

•55 

.l400 

— 

— 

— 

•59 

.1594 

— 

— 

— 

•63 

•1799 

— 

— 
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TABLE I 


- GEOMETRY AKD PRESSURE DATA FOR CIRCULAR CYLIRDERS - Continued 

(c) 7 = 5/3 


M = 00 1 

II 

w 

pq 

Rb/Rs = O.4156J A/Rg = 0 . 2428 ; pg^ = 

0.8805 (0.88132) 

Shock 

Body 

Pi*essure 

ratio 

Sonic 

line 

y 

X 

X 

90° - 

P/Pst 

X 

y 

0 

0 

0,2428 

0 

1.000 

0.123 

0.494 

.03 

.0005 

• 2439 

4 .l 4 

.995 

.16 

•503 

.07 

.0025 

.2487 

9.69 

.969 

.20 

.504 

.11 

.0061 

.2576 

15-3 

.924 

.24 

• 493 

•15 

.0113 

.2707 

21.2 

.860 

.28 

.456 

.19 

.0181 

.2887 

27-2 

.776 

.32 

.348 

-23 

.0265 

.3123 

33.6 

.674 

.34 

.267 

.27 

.0365 

.3430 

40.5 

-553 

— 

— 

-31 

.0482 

• 3840 

48.2 

.411 

— 

— 

•35 

. 06 l 4 - 

— 

— 

— 

— 

— 

.39 

.0763 

— 

— 





— 

.^3 

.0929 

— 

— 

— 

— 

— 

.47 

.1111 

— 

— 


— 

— 

_ 

.1309 

— 


— 

— 

— 

M = 20 

Bg = 0 . 08 j Ri,/Rg = 

0 . 4 ll 6 ; 

A/Rg = 0 . 2435 ; Pst 

= 0.8812 (0.88198) 

0 ^ 

0 

.2435 

0 

1.000 

.124 

.496 

.03 

.0005 

.2446 

4.18 

•994 

.16 

•505 

.07 

.0025 

.2495 

9-79 

.969 

.20 

.506 

• 11 

.0061 

.2585 

15-5 

.923 

.24 

.494 

.15 

.0113 

.2718 

21.4 

.858 

.28 

.438 

.19 

.0181 

.2899 

27.5 

.774 

.32 

•352 

.23 

.0265 

.3138 

34.0 

.611 

.34 

.272 

.27 

.0365 

.3450 

4 i.o 

.550 

— 


.31 

. 048 i 

.3864 

48.9 

.409 

— 

— 

.35 

• 06 i 4 

— 

— 

— 

— 

— 

.39 

.0763 

— 

— 

— 

— 

— 

A 3 

.0928 

— 

— 

— 

— 

— 

• 47 

.1109 

— 

— 

— 

— 

— 

•51 

.1307 

— 

— 

— 

— 

— 

M = 10 

Bg = 0.05 


0 . 4029 ; 

A/Rg = 0 . 2460 ; Pst = 

0.8833 (0.88397) 

0 

0 

.2460 

0 

1.000 

.127 

.503 

.03 

.0005 

.2471 

4.27 

.994 

.16 

.510 

•07 

.0025 

.2521 

10.0 

.968 

.20 

-510 

.11 

.0061 

.2613 

15.8 

.921 

.24 

•497 

•15 

.0113 

.2749 

21.9 

.854 

.28 

.460 

•19 

.0181 

• 2935 

28,1 

•767 

.32 

•355 

.23 

.0265 

.3181 

34.8 

.662 

.34 

.272 

.27 

.0365 

•3503 

42.1 

.539 

— 

— 

.31 

. 048 i 

.3433 

50.3 

.396 

— 

— 

.35 

.0613 

— 

— 

— 

— 

— 

.39 

.0762 

— 

— 

— 

— 

— 

.43 

.0927 

— 

— 

— 

— 

— 

.47 

.1108 

— 

--- 

— 

— 

— 

.51 

.1305 

— 

— 


— 

— 




TABLE I.- GEOMETRY MD PRESSURE DATA FOR CIRCULAR CYLINDERS 

(c) 7 = 5/3 - Continued 


Continued 


M = 8 


Bb = 0; Rb/Eg = 0.3936; A/Rs 



0 

■ZhQj 
■2538 
.2632 
.2772 
.2963 
.3216 
• 35‘^7 
.3990 


90° - Qb 
0 

1^.37 

10.2 

16.2 
22.4 
28.9 
35.8 

43.3 

52.0 


; pg-t = 0.8849 (0.88547) 

Pressure 

ratio 

Sonic 

line 

p/Pst 

X 

y 

1.000 

0.130 

0.509 

.994 

.16 

.515 

•967 

.20 

.515 

.917 

.24 

.502 

.850 

.28 

.468 

.762 

.32 

.368 

.656 

.34 


.533 

— 

— 

.393 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

1 — 

— 

— 

L — 

— 


Bg — -0.02; Rb/Rg 


0.3831; A/Rs = 


0.2515; Pst = 0.8885 (0.88871) 



0 

.2515 

0 

1.000 

.134 

.518 

.03 

.0004 

.2526 

4.49 

.994 

.16 

.522 

.07 

.0024 

•2579 

10.5 

.965 

.20 

.518 

.11 

.0060 

.2676 

16.7 

.915 

.24 

.504 

.15 

.0112 

.2820 

23.1 

.844 

.28 

.463 

.19 

.0180 

.3018 

29.7 

.752 

.32 

.349 

.23 

.0264 

.3282 

36.9 

.641 

.34 

.261 

.27 

.0364 

.3630 

44.8 

.511 


— 

.31 

.o48o 

.4103 

54.0 

.362 

— 

— 

.35 

.0612 

— 

— 

— 

— 



.39 

.0760 

— 

— 

— 

— 

— 

.43 

.0924 

— 

— 

— 

--- 


.47 

.1103 

— 

— 

— 

— 

— 

.51 

.1299 



— 

— 

— 


-0.20; Rb/Rs = 0.3455, 
.2608 

.0004 .2621 

.0024 .2680 

.0060 .2787 

.0112 .2949 

.0180 .3175 

.0264 .3484 

.0363 .3912 

.0478 
.0609 
.0755 
.0916 
•1093 

.1284 
.1490 


M = 4 

0.3455; = 



0.2608; Pst ■ 

I 1.000 f 

.992 
.958 
.897 
.810 
.699 
• 563 
.405 


= 0.8970 (0.89811) 
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TABLE I.- GEOMETRY AJSDD PRESSURE DATA FOR CIRCULAR CYLIHDERS 

(c) 7 - 5/3 - Concluded 


Concluded 


M = 3 

Bs = - 0 . 3 j Rb/Bs = O.313IJ VBs = 0.2728; pst = 0.9116 (0.91153) 


Shock 


Body 


Sonic line 


y 

X 

0 

0 

.03 

.0004 

.07 

.0024 

.11 

.0060 

•15 

.0112 

.19 

.0180 

.23 

.0263 

.27 

.0363 

.31 

.0477 

.35 

.0607 

.39 

•0752 

.43 

.0912 

.47 

.1087 

.51 

.1276 

•55 

.l 48 o 

.59 

.1697 


Bg = - 2.0 

11 

• 

0 

0 

.03 

.0004 

.07 

.0024 

.11 

.0060 

•15 

.0111 

.19 

,0177 

.23 

.0258 

.27 

•0352 

.31 

.0459 

.35 

•0579 

.45 

.0927 

• 55 

.1334 

.65 

•1792 

•75 

.2289 

•85 

.2818 

•95 

.3374 

•99 

.3603 


0 

5.49 

12.9 

20.6 

28.6 

35.8 

47.3 


ratio 


p/Pst 

X 

y 

1.000 

0 .l 64 

0.581 

.991 

.20 

.566 

.950 

.24 

•534 

.878 

.28 

.468 

.774 

.30 

• 393 

.637 

.32 

.280 

.465 

.34 

.208 

— 

— 



= 0.1938; /VRs = 0.2849; Pst = 0.9504 (0.9518) 


.2849 

.2872 

•2979 

.3189 


0 

8.90 

21.2 

34.6 


.3564 50.7 


1.000 

.980 

.888 

.716 

.456 



1.00 

.940 

.675 
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(b) w(s,t) as function of t. 


A 


Figure 1.- Concluded 


rn oj 



o)(s. 



(a) Modified stream function w[ (l/2) At] . 

Figure 2.- Effect of smoothing procedure on the calculations 
M = 7 = 1*J+, Bg = -0.06. 











ai(s,t) 



(a) iA)(s^t) vs. t for w[0_, (l/2)At] = 1.1000. 


Figure 3*“* Effect of artificial error on the smoothed calculations 








(c) 0.2l) vs. s for w[0^(l/2)At] = 1.1000 


Flgixre 3'“ Concluded. 












(d.) B]2 vs. Bg • 
Figure 4.- Concluded. 



kk 



X 


Figure 5-“ Effect of parameter A on body shape; M = 10 , 7 = 1.4, Bg = 0 . 1 . 
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